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^ | Abstract 

[t ! | We integrate nonabelian Toda field equations [Kr] for root systems of types A, 

B, C, for functions with values in any associative algebra. The solution is expressed 
via quasideterminants introduced in [GR1],[GR2], [GR4]. In the appendix we re- 
view some results concerning noncommutative versions of other classical integrable 
equations. 

Introduction 

Nonabelian Toda equations for the root system A n -\ were introduced by Polyakov 
(see [Kr]). They are equations with respect to n unknowns (p — (0i> •••> 0«) £ 
! A[[u, v]], where A is some associative (not necessarily commutative) algebra with 
*b5> unit: 

> ■ K J du V dv ^ 

Suppose that A is a *-algebra, i.e. it is equipped with an involutive antiautomor- 
phism * : A —>■ A. Then, setting in (1) n +i-i = (0*) _1 , we obtain a new system 
of equations. If n = 2k, we get the nonabelian Toda system for root system 




d 



>>2<t>i\ 3 = 1 



[2) i- ( = h+tfj 1 - MJ^ 2<3<k-l 



ty- 1 *!; 1 - Mali, 



j = k 



If n = 2k + 1, we get the nonabelian Toda system for root system B^\ 



W 1 , 



(3) #- ( ^fc 1 ) = I - M7\> 2<3<k 



where 



du V 0u J 7 J J 

■ l j,* j i-i 



i = i 

2<J 

A + l - .'A • i o k 1 . j = k + 1, 



Quasideterminants were introduced in [GR1], as follows. Let X be an m x m- 
matrix over A. For any 1 < i,j < to, let r,(X), Cj(X) be the i-th row and the 
j-th column of X. Let X lJ be the submatrix of X obtained by removing the i-th 



row and the j-th column from X. For a row vector r let be r without the j-th 
entry. For a column vector c let be c without the i-th entry. Assume that J y 
is invertible. Then the quasideterminant \X\ij G A is defined by the formula 

(4) \X\ ij = x ij -r i (X)^(X^)- 1 c j (X)( i \ 

where the ij-th entry of X. 

In this paper we will use quasideterminants to integrate nonabelian Toda equa- 
tions for root systems of types A, B, C. We do not know yet how to generalize these 
results to root systems of types D-G, and to affine root systems. 

Our method of integration, which is based on interpreting the Toda flow as a flow 
on the space of factorizations of a fixed ordinary differential operator. This method 
and explicit solutions of Toda equations are well known in the commutative case 
(see [LS]; see also [FF] and references therein; for the one variable case see [Ko]). 
1. Factorization of differential operators 

Let R be an associative algebra over a field k of characteristic zero, and D : R — > 
R be a fc-linear derivation. 

Let fi, f m be elements of R. By definition, the Wronski matrix W(fi, f m ) 

is 

/ fl ••• fm \ 

Dh - Df m 



W(f 1 ,...J m ) = 



KD™- 1 ^ ... D rn ~ 1 f m J 



n > 



We call a set of elements /i, f m eR nondegenerate if W(fi, f m ) is invertible 

Denote by R[D] the space of polynomials of the form aoD n + ai-D n_1 + ... + a 
cii G R. It is clear that any element of R[D] defines a linear operator on R. 

Example: R = C°°(]R), D = ^. In this case, R[D] is the algebra of differential 
operators on the line. 

By analogy with this example, we will call elements of R[D] differential operators. 

We will consider operators of the form L = D n + aiD n_1 + ... + a n . We will call 
such L an operator of order n with highest coefficient 1. Denote the space of all 
such operators by R n {D). 

Theorem 1.1. (i) Let fi,...,f n G i? 6e a nondegenerate set of elements. Then 
there exists a unique differential operator L G R[D] of order n with highest coeffi- 
cient 1, such that Lfi = for i = 1, ...,n. R is given by the formula 

(1-1) Lf=\W(f 1 ,...J n J)\ n+1 , n+1 . 

(ii) Let L be of order n with highest coefficient 1, and /i,...,/ n be a set of 
solutions of the equation Lf = 0, such that for any m < n the set of elements 
/i,...,/ m is nondegenerate. Then L admits a factorization L — (D — b n )...(D — b\), 
where 

(1.2) b i = (DW i )W~ 1 , W i = \W(f 1 ,...,f i )\ ii . 

Proof, (i) We look for L in the form L = D n + aiD n ~ 1 + ... + a n . From the equations 
Lfi = it follows that 



(a n , 0l ) = -(D n f u D n f n )W{h, f n )-\ 



By definition, 

\W(f U -, fn, f)\n+l,n+l = D n f - (D n f lt D n f n )W(h, f n )~\f, Df, D n ~ 1 f) T = 

D n f + (a n , oi)(/, Df, D n - x f) T = Lf. 

(ii) We will prove the statement by induction in n. For n = 1, the statement is 
obvious. Suppose it is valid for the differential operator L n _i of order n — 1 with 
highest coefficient 1, which annihilates /i, / n _i (by (i), it exists and is unique). 
Set b n = (DW n )W~ l , and consider the operator L = (D — b n )L n _i. It is obvious 
that Lfi = for i = 1, n — 1. Also, by (i) 

L/ n = (D - b n )L n _ x f n = (D- b n )W n = 0. 

Therefore, by (i), L = L. □ 

Now consider the special case: R = A[[t]], where A is an associative algebra over 
k : and D = A (here t commutes with everything). In this case, it is easy to show 
that nondegenerate sets of solutions of Lf = exist, and are in 1-1 correspondence 
with elements of the group GL n (A), via f = (fi, f n ) — > VF(f)(0). 

It is clear that two different sets of solutions of the equation Lf = can define 
the same factorization of L. However, to each factorization 7 of L we can assign 
a set f 7 = (/1, f n ) of solutions of Lf = 0, which gives back the factorization 7 
under the correspondence of Theorem l.l(ii). This set is uniquely defined by the 
condition that the matrix W(fj)(0) is lower triangular with 1-s on the diagonal. 

Here is a formula for computing f 7 , which is well known in the commutative 
case. 

Proposition 1.2. If 7 has the form 

L = (D- (Dg^g- 1 )...^ - (D 9l )g^), 
where gi(0) = 1, then f 7 = f n ), where 

(1.3) fj(t)= [ I' ... I' 2 g 1 (t)g 1 (t 1 )- 1 g 2 (t 1 )g 2 (t 2 )- 1 ...g J (t J - 1 )dt J . 1 ...dt 2 dt 1 , 
Jo Jo Jo 

where J' Q U (Y1 a % t l )dt := J2 a ii^T ■ 

Proof. It is easy to see that if f = / n ), with fj given by (1.3), then W(f)(0) is 

strictly lower triangular. So it remains to show that fj is a solution of the equation 
Ljf = 0, where L, = (D — [ I)g l )g J 1 )...[ I) - (D 9l )g^). 

We prove this by induction in j. The base of induction is clear, since from (1.3) 
we get fi = g\. Let us perform the induction step. By the induction assumption, 
From (1.3), we have 

/;(*) = Si (*) / 9\{s)- x h{s)ds, 
Jo 

where h obeys the equation (D — (Dgj)g~ 1 )...(D — (Dg 2 )g 2 ~ 1 )h = 0. Thus, we get 



(D-(Dg 1 )g- l )f J =h. 



This proves that Ljfj = 0. □ 

Now consider an application of these results to the noncommutative Vieta the- 
orem [GR3]. Let A be an associative algebra. Call a set of elements 
generic if their Vandermonde matrix V(a;i, ...,x n ) (V^ := x l ~ ) is invertible. 

Consider an algebraic equation 

(1.4) x n + aix n - 1 + ... + a n = 0. 

with aj G A. Let xi,...,x n G A be solutions of (1.4) such that x±,...,Xi form a 
generic set for each i. Let V(i) = V( x i> ...,Xi), and y t = | V^(i) | ^0;^ j V (i) ] ^ . 

Theorem 1.3. [GR3] (Noncommutative Vieta theorem) 
(1-5) a r = (-l) r Vir-Vi!- 

i 1 <...<i r 

Proof. (Using differential equations.) Consider the differential operator with con- 
stant coefficients in R = A[[t]]: 

L = D n + a 1 D n ~ 1 + ... + a n . 

We have solutions /j = e tXl of the equation Lf = 0, and for any i the set /i, /j 
is nondegenerate, since its Wronski matrix W(i) is of the form 

(1.6) W(») = V(i)diag(e tx \...,e tXi ). 

Thus, by Theorem 1(h), the operator L admits a factorization 

L=(D-b n )...(D-h), 

where hi = (D\W {i)\u)\W {i)^ . Substituting (1.6) into this equation, we get 

(1-7) b i = \V{i)\ ii x i \V{i)\7 1 = y i . 

Since [D, 6j] = 0, we obtain the theorem. 

2. Toda equations 

In this section we will solve equations (1),(2),(3). 

Let M = {{g,h) G A[[t]] © A[[t]] : g(0) = h(0) G A*}, where A* is the set of 
invertible elements of A. We have the following obvious proposition, which says 
that solutions of Toda equations are uniquely determined by initial conditions. 

Proposition 2.1. The assignment <fi(u,v) — > (4>(u, 0), 0(0, v)) is a bijection be- 
tween the set of solutions of the Toda equations (1) and the set M. 

Let B be an algebra over k, R = B[[v]], D = -J^. For <p = ...,0 n ), where 
4>i G B[[v]] are invertible, define Lf G R[D] by 

(2.1) Lf = (L> - (D(f) i )(f)~ 1 )...(D - (Dcp^- 1 ). 

Now set B = A[[u]]. Then = A[[u, v]], so for any (j) = (0i, ...,0 n ), with all 

0, G A[[it, v]] invertible, we can define Lf . 



Proposition 2.2. A vector-function (f>(u,v) is a solution of Toda equations (1) if 
and only if 

(2-2) ^ = -0 l+1 0- 1 Lf_ 1 ,z<n-l; = 0. 

Proof. Let be a solution of the Toda equations. Set bi = {D$i)<f>~ x . We have 
Lf = Li = (D — bi)...{D — b\). Therefore, for i < n — 1 we have 

8T 1 

= - J2(D - h)...{D - bj+Jh+KfiJ^D - &,_!)...(£> - h) 

3 = 1 
i-1 

(2.3) + J2(D ~ bi)...(D - bj+^+KPj^D - bj)...(D - h). 

i=i 

But 

(D - b j+1 ) o (Pj+kPJ 1 = 4> j+ i^-\D - b 3 ). 

Therefore, the right hand side of (2.3) equals — <f>i + i(j)~ Li_\. The same argument 
shows that for i = n the derivative ^j 2 - vanishes. 

Conversely, it is easy to see that equations (2.2) imply (1). □ 

Now we will compute the solutions of Toda equations explicitly. Let 771, rj n , ipi, ip n G 
A[[t}} be such that 7^(0) = ipi(0) G A*. We will find the solution of the following 
initial value problem for Toda equations: 

(2.4) fofa 0) = Vi(u), 0i(O, v) = rn(v) 

Proposition 2.1 states that this solution exists and is unique. 

Let gi(v) = ?7i(i')^i(0) _1 . Let f = (/1, ...,f n ) be given by formula (1.3) in terms 
of gi. Define the lower triangular matrix A (it) whose entries are given by the 
formula 

Ay(ti) = / / ... / ^(ti-j^i^U- j)^i-i(ti-j -1)... ip~ 1 (t 1 )ip j (u)dt i - j ...dt 1 

Jo Jo Jo 

Let P = (fi, f%) be defined by the formula 

P = fA(«). 

Then we have 

Theorem 2.3. The solution of the problem (2.3) is given by the formula 
(2-5) cf )i (u,v) = \W(f?,...J?)\ii 

Proof. Let L" = L^ u,v \ where (f)(u,v) is defined by (2.5). By Theorem l.l(ii), we 
have L^fJ = for j < i. Differentiating this equation with respect to u, we get 



(2.6) 



On the other hand, it is easy to see that 



where 



This implies that 



A'(u) = A(u)Q(u), 

(V^V)^O) i = j 
1 i=j+l 
otherwise 



QfU 



df u 

Therefore, L^-^- = for i > j + 1, and 



Thus, from (2.6) we get 



By Theorem l.l(i), this implies equations (2.2), which are equivalent to Toda equa- 
tions. 

It is obvious that 4>(u, v) satisfies the required initial conditions. The theorem is 
proved. □ 

If initial conditions (2.4) satisfy the symmetry property (j) n+ i-i = (0*) _1 , then 
we obtain a solution of the initial value problem for equations (2) (for even n), 
and (3) (for odd n). This gives a complete description of solutions of systems of 
equations (1),(2),(3). 

In the case (p(u, 0) = 1, the Toda flow can be interpreted as a flow on the space 
of factorizations of a differential operator. Namely, let L be a differential operator 
of order n with highest coefficient 1. Let F(L) be the space of factorizations of L. 
Let N~ be the group of strictly lower triangular matrices over A. It is easy to see 
that the map it : F(L) — > N~ , given by = H / (f T )(0), is a bijection. We will 
identify F(L) with N~ using 7r. 

Let j(u) = j(0)e uJn be a curve on N~ generated by the 1-parameter subgroup 
e" J ™, where J n is the lower triangular nilpotent Jordan matrix ((J n )ij = ^j+i)- 
Let L = (D — b™)...(D — b\) be the factorization of L corresponding to the point 
7(it). Let (pi(u,v) be such that (Dcpi)^ 1 = 6", and (pi(u,0) = 1. (here D = ^). 
Then we have the following Corollary from Theorem 2.3. 

Corollary 2.4. </> = (0 l5 ...,0 n ) is a solution of Toda equations (1), and all solu- 
tions with 4>(u, 0) = 1 are obtained in this way. 

Proof. For the proof it is enough to observe that if ipi(u) = 1 then Q(u) = J n , and 
A(u) = e uJn . □ 

Remark. An analogous statement can be made for equations (2) and (3). In this 
case, instead of an arbitrary differential operator L e R n (D) one should consider 
a selfadjoint (respectively, skew-adjoint) operators, instead of the group N~ - a 



maximal nilpotent subgroup in Sp n (A) (respectively, n (A)), and instead of J n the 
sum of simple root elements in the Lie algebra of this group. In the commutative 
case, such a picture of the Toda flow is known for all simple Lie groups [FF]. 
Consider now the statement of Theorem 2.3 for % = 1. In this case we have 

n 

(2-7) Mu,v) = A» = J2fi(v)\i(u) 

i=i 

Thus, we get 

Corollary 2.4. [RS] If (f) is a solution of the Toda equations then (f)\(u,v) = 
Yl7=i Pi( u ) ( li( v )> where pi, qi are some formal series. 

Now let us discuss infinite Toda equations. These are equations (1) with n = oo. 
These equations allow to express (pi recursively in terms of / = (pi, which can be 
done using quasideterminants: 

(2.8) j> i = \Y i (f)\ ii , Yi {f):=w(j,^,...,^y 

where are regarded as functions of v for a fixed u. This formula is easily proved 
by induction. It appears in [GR2], section 4.5 in the case u = v, and in [RS] in the 
2-variable case. 

Formula (2.8) can be used to give another expression for the general solution of 
finite Toda equations, which appears in [RS]. Indeed, we have the following easy 
proposition. 

Proposition 2.5. Let f e A[[u, v]] be such that Y 1 (f), ...,Y n (f) are invertible ma- 
trices. Then |Y^ + i(/)| n+ i n+ i = if and only if f is "kernel of rank n", i.e. 

n 

(2-9) / = X>* 

i=i 

Thus, taking / = <pi of the form (2.9), with Yi,...,Y n invertible, and using 
formula (2.6) we will get a solution (0 1; n ) of the finite Toda system of length 
n. It is not difficult to show that in this way one gets all possible solutions. 

Example. Consider the nonabelian Liouville equation 

which is a special case of (2) for k = 1. Consider the initial value problem 

(2.11) 0(0, v) = v(v),<i>(u, 0) = ip(u),r](Q) = ^(0) = a. 
By Theorem 2.3, we get the following formula for the solution: 

(2.12) 4>(u,v) =n(v)(a- 1 +J^ J r ] - 1 {t){r ] *)-\t)a%^)- 1 {s)^-\s)dsd^{u). 



For example, if i/j(u) = 1, we get 



(2.13) <f>(u,v)=ri(v)(l + uj rj-WtfrWdty 

In the commutative case, these formulas coincide with the standard formulas for 
solutions of the Liouville equation. 

Appendix: Noncommutative soliton equations 

In this appendix we will review some results about noncommutative versions of 
classical soliton hierarchies (KdV,KP). These results are mostly known or can be 
obtained by a trivial generalization of the corresponding commutative results, but 
they have never been exposed systematically. 

We will follow Dickey's book [D]. 

Al. Nonabelian KP hierarchy. 

Nonabelian KP hierarchy is defined in the same way as the usual KP hierarchy 
[D]. Let 

(Al) L = d + wod' 1 + Wl d~ 2 + .... 

be a formal pseudodifferential operator. Here d = Wi G -A [[#]], where A is an 
associative (not necessarily commutative) algebra with 1. Consider the following 
infinite system of differential equations: 

rlT 

(A2) — - = [B m ,L],B m = (L m )+, 

Otm 

where for a pseudodifferential operator M, we denote by M + ,M_ the differential 
and the integral parts of M (i.e. all terms with nonnegative, respectively negative, 
powers of d). For brevity we will write L± instead of (L m )±. 

Each of the differential equations (A2) defines a formal flow on the space P of 
pseudodifferential operators of the form (Al). Indeed, [L™, L] = — [L™, L], and the 
order of [L™, L] is at most zero. 

Proposition Al. The flows defined by equations (A2) commute with each other. 

Proof. The proof is the same as the proof of Proposition (5.2.3) in [D]. Namely, one 
needs to check the zero curvature condition 

< A3 > ff-f 2 -^^ - 



which is done by a direct calculation given in [D]. 

The hierarchy of flows defined by (A2) for m = 1,2,3, ... is called the noncom- 
mutative KP hierarchy. 

A2. Noncommutative KdV and nKdV hierarchies. 

As in the commutative case, we can restrict the KP hierarchy to the subspace 
P n C P of all operators L such that L n is a differential operator, i.e. = [D]. 

This subspace is invariant under the KP flows, as -q^- = [L™, L"]_. The space P n 
can be identified with the space D n of differential operators of the form 

(A4) M = d n + u 2 d n - 2 + ... + u n , 



by the map L — > M = L n . The KP hierarchy induces a hierarchy of flows on D n 
called the nKdV hierarchy: 

(A5) ^ = [M? /n ,M}. 

Among these flows, the flows corresponding to m = nl, I e Z + , are trivial, but the 
other flows are nontrivial. 

For n = 2 the nKdV hierarchy is the usual KdV hierarchy. The first two non- 
trivial equations of this hierarchy are 

(A6) u tl = u x , u ta = |(« ra x + 3u x u + 3uu x ). 

The second equation is the noncommutative KdV equation. 
A3. Finite-zone solutions of the nKdV hierarchy. 

Denote the vector fields of the noncommutative nKdV hierarchy by V m , m ^ nl. 
Let V = Ym=i a iVi> a p 7^ 0, be a finite linear combination of Vi with coefficients 
Oi G I. Let S(a), a = (oi,...,o p ) be the set of stationary points of V in D n = 
Affx]]"- 1 . 

The set S(a) can be identified with A^ -1 ^. Indeed, S(a) is the subset of D n 
defined by the differential equations V{M) = 0, which have the form 

(A7) f£ =*!(«,,«},...,,,?-'>). 

So Ui are uniquely determined by their first p derivatives at 0. 

Since the vector field V commutes with the nKdV hierarchy, S(a) is invariant 
under the flows of this hierarchy. Using the identification S(a) — > A^ n ~ 1 ^ p , we can 
rewrite each of the nKdV flows as a system of ordinary differential equations on 
A^ n ~ r>p . Thus, solutions of nKdV equations belonging to S(a) can be computed 
by solving ordinary differential equations. Such solutions are called finite-zone 
solutions. 

Remark. If M e S(a) then [Q(M 1 / n ) + ,M] = 0, where Q(x) = J2^x\ Thus, 
we have two commuting differential operators M and = Q(M 1 ^ n ) + in 1 variable. 
If the algebra A is finite dimensional over its center then there exists a nonzero 
polynomial R(x, y) with coefficients in the center of A, such that R(M, N) = 0. 
This polynomial defines an algebraic curve, called the spectral curve. The operator 
M can then be computed explicitly using the method of Krichever [Kr2] . Similarly, 
all nKdV equations resticted to the space of such operators can be solved explicitly 
in quadratures. However, if A is infinite-dimensional over its center, the polynomial 
R does not necessarily exist, and we do not know any way of computing M explicitly. 

A4. Multisoliton solutions of the noncommutative KP hierarchy. 

Here we will construct N-soliton solutions of the KP and KdV hierarchies in the 
noncommutative case. We will use the dressing method. In the exposition we will 
closely follow Dickey's book [D]. 

We will now contruct a solution of the KP hierarchy. Let t = (ti,t 2 , ...) Consider 
the formal series 



(A8) 



£(cc, t, a) = (x + t\)a + t20? + ■■■ + t r a r + a E A. 



Fix cki, ocn, pi, (3n, a>i, £ A and set 

(A9) y fl (M) = e^'*' 01 ^ + a a e^ x '*' p '\ 

Define the differential operator of order iV, with highest coefficient 1, by 

(A10) = |W(yi,...,yAr,/)|iv+i^+i(a;), 

where the derivatives in the Wronski matrix are taken with respect to x, and we 
assume that the functions y\, j/at are a generic set (in the sense of Chapter 1). 
Set 

(All) L = $d$ _1 . 

Proposition A2. The operator-valued series L(x, t) is a solution of the KP hier- 
archy. 

Proof. The proof is the same as the proof of Proposition 5.3.6 in [D]. 
Such solutions are called N-soliton solutions. 

Proposition A2 can be used to construct N-soliton solutions of the nKdV hierar- 
chy. Namely, we should restrict the above construction to the case when (3k = £k&k, 
where is an n-th root of unity. In this case it can be shown as in [D] that the 
operator M = L n = $9 n $ _1 is a differential operator of order n, and it is a solution 
of the nKdV hierarchy. 

As an example, let us consider the N-soliton solutions of the KdV hierarchy. In 
this case, we may set t^k = 0, and we have 

(A12) y s = e^ x ^ aa) + a s e~^ x ^ aa) . 

Proposition A3. Let 

(A13) b i = (dW i )W~ 1 , W i :=\W(y 1 ,...,y i )\ ii . 

Then the function 

N 

(A14a) u(x,t) = 2d(^bi). 

i=i 

is a solution of the noncommutative KdV hierarchy. 

Proof. Let $ = d N + vid N ~ x + .... Then from the equation (<9 2 + u)$> = we 
obtain u = —2dv\, and from Theorem 1.1(h) v\ = — J^&i, where hi are given by 
(A13). This implies (A14a). 

Another formula for u, which is equivalent to (Al4a), is the following. Let 
Y(yi, yjv) be the matrix which coincides with the Wronski matrix W(yi, yjv) 
except at the last row, where it has y^ instead of y\ N . Let Yn = \Y(yi, yN)\NN- 
Then the function u given by (Al4a) can be written as 



(Al4b.) 



u = 2d(Y N W^), 



In the commutative case formulas (Al4a),(Al4b) reduce to the classical formula 



(A15) u = 2d 2 \n[fetW{ yi ,...,y N )}. 

In particular, we can obtain N-soliton solutions of the noncommutative KdV 
equation u t = j(u xxx + 3u x u + 3uu x ). For this purpose set ti = 0, i 7^ 3, and 
£3 = £. Then we get 

(A16) y s = e QsX+Q ^ + a 8 e- asX -°^\ 

and the solution u(x,t) is given by (Al4a),(Al4b). 

For example, consider the 1-soliton solution. According to (Al4a), it has the 
form 

(A17) u = 2^-[{e ax+aH - ae- ax - aH )a(e ax+aH + ae- ax - aH )- 1 ]. 

ox 

In the commutative case, it reduces to the well known solution 

- 2C " 2 - 1 In 

cosh 2 (ax + aH — c) ' 2 

- the solution corresponding to the solitary wave which was observed by J.S.Russell 
in August of 1834. 
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